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ABSTRACT 

Runway growth is an important stage in planet formation during which large protoplanets form, 
while most of the initial mass remains in small planetesimals. The amount of mass converted into 
large protoplanets and their resulting size distribution are not well understood. Here, we use analytic 
work, that we confirm by coagulation simulations, to describe runaway growth and the corresponding 
evolution of the velocity dispersion. We find that runaway growth proceeds as follows: Initially all the 
mass resides in small planetesimals, with mass surface density a, and large protoplanets start to form 
by accreting small planetesimals. This growth continues until growth by merging large protoplanets 
becomes comparable to growth by planetesimal accretion. This condition sets in when E/cr ~ a'^/^ ^ 
10^"^, where E is the mass surface density in protoplanets in a given logarithmic mass interval and a is 
the ratio of the size of a body to its Hill radius. From then on, protoplanetary growth and the evolution 
of the velocity dispersion become self-similar and E remains roughly constant, since an increase in E by 
accretion of small planetesimals is balanced by a decrease due to merging with large protoplanets. We 
show that this growth leads to a protoplanet size distribution given by N{> R) oc where N{> R) 
is the number of objects with radii greater than R (i.e., a differential power-law index of 4). Since 
only the largest bodies grow significantly during runaway growth, E and thereby the size distribution 
is preserved. We apply our results to the Kuiper Belt, which is a relic of runaway growth where planet 
formation never proceeded to completion. Our results successfully match the observed Kuiper belt 
size distribution, they illuminate the physical processes that shaped it and explain the total mass that 
is present in large Kuiper belt objects (KBOs) today. This work suggests that the current mass in 
large KBOs is primordial and that it has not been significantly depleted. We also predict a maximum 
mass-ratio of Kuiper belt binaries that formed by dynamical processes of a~^/^ ~ 10, which explains 
the observed clustering in binary companion sizes that is seen in the cold classical belt. Finally, our 
results also apply to growth in debris disks, as long as frequent planetesimal-planetesimal collisions 
are not important during the growth. 

Subject headings: circumstellar matter — Kuiper belt — planetary systems: general — planets and 
satellites: formation — solar system: formation 



1. INTRODUCTION 

Ice giants, the cores of gas giants and protoplanets 
that later form terrestial planets, are generally believed 
to have formed by coagulation from small planetesi- 
mals. Understanding the evolution of the size distri- 
bution of growing protoplanets, their velocity dispersion 
and the interplay between the two is crucial for shed- 
ding light onto the planet formation process. Here, we 
study the runaway growth of protoplanets (e.g . Safronov 
1973 IGreenberg et al.l [l978l: IWetherih fc Stewart .1989. 



Kokubo fc Idalll996[ l and their subsequent velocity evo- 
lution. Runaway growth can occur when the accretion 
cross section of protoplanets is enhanced by gravitational 
focusing. Since, gravitational focusing is strongest for the 
largest bodies, the radii of larger protoplanets run away 
from that of smaller. This leads to a size distribution that 
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develops a tail containing a small number of large proto- 
planets. It is however not always clear what fraction of 
the total mass participates in this runaway growth and 
what protopla net size distributi on such ru naway growth 
gives rise to (|Lei [2000: Malvs hkin fc Goo dman 2001). 
In this paper we address both of these questions. In 
the following sections, we discuss runaway growth in the 
context of the Kuiper belt, which is an ideal laboratory 
to test our results, since it is a remnant of the primor- 
dial Solar system, where planet formation never reached 
completion. The results however, also apply to runaway 
growth during planet formation and growth in debris 
disk, as long as gas plays no significant role in the ac- 
cretion and damping of the velocity dispersion. 

The Kuiper belt consists of a disk of icy bodies lo- 
cated at the outskirts of our planetary system, just 
beyond the orbit of Neptune and contains some of 
the least processed bodies in our Solar system. Mo- 
tivat ed by the discovery of the first Kuiper belt ob- 
ject (jJewitt fc Luul [19931 ) after Pluto and Charon, sev- 
eral groups conducted large scale surveys to character- 
ize the Kuiper belt. These efforts led to the discov- 
ery of more than 1200 objects in the Kuiper belt to 
date. The Kuiper belt size distribution contains many 
important clues concerning the formation of Kuiper belt 
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objects (KBOs), their effective s t rength and their colli- 
sional evolution (lpohiianvilll969HStern fc Colwelil 119971: 
Davis fc Farinenall997HKenYon fc Luulll999l:lPan fc Saril 
20051 ). It also provides a snapshot of an earlier stage 
of planet formation, which was erased elsewhere in the 
Solar system where planet formation proceeded all the 
way to completion. The cumulative size distribution of 
KBOs larger than i? > 50 km is well described by a single 
power- law given by 

N{>R)(xR^-'' (1) 

where A^(> R) is the number of objects with radii greater 
than R, and q is the power-law index. Kuiper belt 
surveys find that the size distribution for KBOs with 
radii greater than about 50 km follow s this power-law 
with Q ^ A (e.g. 'Truiillo ct al." 2001'; 'Bernstei n et al.l 
l2004HFuentes fc Holm an 2008; Fraser.et _al. 2008), which 
implies roughly equal mass per logarithmic mass in- 
terval. This size distribution is a relic of the ac- 
cretion history in the Kuiper belt and therefore pro- 
vides valuable insights into the formation of large KBOs 
(R > 50 km) (e.g. Stern 1996; Davis & Farinclla 1997; 
iKenyonI 120021 ) . Observations suggest that there is a 
break in the powe r-law size distributio n at s maller KB O 
sizes (e.g. [B ernstei n et "aIll2004l;lFuente s fc H olmanll2008l; 
Fraser et al.li2008l; iSchUchting et aLn2009t iFuentes et al l 
20101 ). The break in the size distribution is generally at- 
tributed to collisions that break-up small KBO s (i.e.,/? < 
50 km ) and modify their size di s tribution (e. g. IDohnanvil 
[19691: IKenvon fc Bromlevll200l Wa^k Saril I2005D . The 
KBO size distribution below the break is still poorly 
constrained, although some encouraging progress has 
been made recently in probing the abund ance of sub-km- 
sized KBOs by s tellar occultations (e.g. | Liu et al. 2008; 
iSchlichting etaIll2009HBianco et al.il20m 

The work presented in this paper focuses on the size 
distribution of large KBOs (/? > 50 km), which is well 
constrained by observations and which sheds light onto 
the formation of KBOs, protoplanets and accretion pro- 
cesses that could be ongoing in other debris disks. Nu- 
merical coagulation simulations have been successful in 
reproducing the observed KBO size distribution. Such 
simulations typically find that the accretion processes 
of KBOs yield a power-law size distrib ution with q ~ 
3.8—4.5 for 10-100 km and larger objects (Ken von fc Lud 
II999; Kenyon 2002; Kcnyon & Bromley 2004), which is 
consistent with the observed power-law size distribution. 
Despite their success, the reason for the actual slope of 
the distribution has so far not been explained by such 
simulations. In this paper, we offer an explanation for 
the slope of the KBO size distribution and for the amount 
of mass in the large KBOs that are observed in today's 
Kuiper belt. Specifically we find a power-law index of 
g ~ 4 and a total mass in large KBOs of ~ lO^'^ of 
the initial planetesimal mass, which is consistent with 
the current observed mass in the Kuiper belt. We also 
make a prediction for the maximum mass ratio of Kuiper 
belt binaries that formed by dynamical processes, for ex- 
ample, by three body interactions, and show that our 
prediction is in good agreement with the observations. 
Although our work focuses on the Kuiper belt, the re- 
sults also apply to early stages of planet formation and 
protoplanetary growth in debris disks. 



Our paper is structured as follows: In §2 we analyti- 
cally describe the growth of large KBOs, including their 
velocity dispersion and derive the slope of the KBO size 
distribution. We confirm our analytic results in §3 with 
coagulation simulations. In §4, we discuss how semi- 
collisional accretion, binary mergers and frequent plan- 
etesimal collisions would affect our results. We show that 
our results on the KBO growth and velocity dispersion 
have interesting implications for the formation of Kuiper 
belt binaries and predict the maximum mass ratio for 
binaries that formed by dynamical processes in §5. Dis- 
cussion and conclusions follow in §6. 

2. ANALYTICAL TREATMENT 

In order to gain an analytic understanding of the 
growth processes of large KBOs and the associated ve- 
locity evolution we use the 'two-groups approximation' 
()Goldreich et al . 2002, 2004b). The 'two-groups approx- 
imation' consists of the identification of two groups of 
objects, small ones, that contain most of the total mass 
with mass surface density ct, and large ones, that contain 
only a small fraction of the total mass with mass surface 
density S ^ ct. We define S as the mass surface density 
in a single logarithmic mass interval, that includes the 
largest bodies formed at a given time. In contrast ct is 
defined as the total mass in small objects. 

Within the framework of the 'two-groups approxima- 
tion' we arrive at the following picture for KBO growth. 
Initially all the mass is in small bodies. As the small 
bodies start to accrete each other, large bodies begin to 
form. To simplify the argument, we only consider the 
mass surface density of the small and large bodies here, 
ignoring intermediate size bodies for now. As we show 
later, the large and small bodies alone determine the ve- 
locity dispersion for bodies of all sizes and only the large 
bodies grow significantly. In the initial stage, S grows 
due to the accretion of small bodies. Therefore the size 
of the largest bodies and the total mass in large bod- 
ies increases with time. During this growth phase the 
velocity dispersion of the small bodies increases due to 
viscous stirring by the large bodies. The velocity dis- 
persion of large bodies is damped by dynamical friction 
provided by the small bodies. E continues to grow un- 
til the growth of large KBOs by accretion of comparable 
size objects starts to compete with growth by accretion 
of small bodies. From then on, E remains roughly con- 
stant in a given logarithmic mass interval, while the size 
of the large KBOs grows linearly with time. How the 
KBO growth ended and how exactly the small bodies 
were lost from the Kuiper belt are still the subject of 
ongoing research and are unimportant for the purpose of 
this paper and we therefore will not discuss them here 
further. We confirm the outlined KBO growth analyti- 
cally and with numerical simulations. We show that the 
mass-ratio, E/ct, is not arbitrary but an outcome of KBO 
growth. 

2.1. Growth and Velocity evolution 

Large KBOs viscously stir the small bodies, increasing 
the small bodies' velocity dispersion u. As a result u 
grows on the same timescale as R, as long as the small 
bodies experience no significant damping by either gas or 
mutual collisions, which are, most likely, not yet impor- 
tant (see section|4]). We can therefore write the evolution 
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of the small bodies' velocity dispersion as 

-4 



1 du 
u dt 



pR \VH 



(2) 



where Q is the orbital frequency around the Sun, R the 
radius of the the large bodies and p is their material 
density. The Hill velocity of the large bodies is vh , which 
is given by vh — VLRh where Rh is the Hill radius, 
Rh = a(M/3MQ)i/3, where Mq is the mass of the Sun, 
M and a are the mass and semi-major-axis of the large 
KBOs. Finally, a — R/Rh and is ~ 10~^ at the distance 
of the Kuiper belt. We assumed in writing the expression 
for u that u > vh, we will verify that this is the correct 
and self-consistent velocity regime for u at the end of this 
section. Initially, the large bodies grow by accreting the 
small ones. Their growth rate is given by 



I dR ^ cr _ 

— ^ \l a 

R dt pR \vh 



(3) 



where we assumed that v < vh, which we verify later in 
this section. Equating the rates from equations ([2]) and 
and solving for u we find 



VH 



1/2 



(4) 



(jGoldreich et al.ll2004bD . The velocity v of large KBOs 
increases due to mutual viscous stirring, but is damped 
by dynamical friction from the sea of small bodies such 
that V < u. The competition between the stirring and 
damping can be written as 



1 dv 
V dt 



pR 



VH 



pR \VH 



(5) 



Balancing the stirring and damping rates for v and sub- 
stituting for u from equation Q into equation ([5]), we 
find 

3 

VH \(J J 

Having derived expressions for the large and small bod- 
ies' velocity dispersions we now turn to examining the 
growth in more detail. 

Expressions similar to the ones a bove have been 
widely u sed in the literature (e.g. ISa fronovl 197^ 
'Grcenber g et al.lll991 : IDones fc Tremaine 1993; Rafikov 
2003t iGoldreich et al.[[2004bt lChiang'eni.ii2007i ). How- 
ever, in these works E/cr was treated as a free parameter. 
Here we present our method for deriving E/cr and calcu- 
late its value during runaway growth. 

Within our 'two-groups approximation' large bodies 
have two distinct modes for growth. In the first, they 
can grow by the accretion of small bodies. In the sec- 
ond, they grow by accreting objects comparable to their 
own size. The growth rate for large KBOs is given by 



1 dR 
R~dt 



u 

VH 



E 

n a 

pR 



3/2 



(7) 



where we used the accretion rate corresponding to sub- 
Hill velocity dispersions for large bodies and the corre- 
sponding rate for u > vh for small bodies. The first term 



in expression ([7]) describes the growth of large KBOs by 
the accretion of small bodies and the second corresponds 
to the growth of large KBOs by accreting objects of their 
own size. Substituting the expression for u from equation 
Q we have 



1 dR 
'R~dt 



E 

'Jr^ 



(8) 



Comparing the two terms in expression ([5]) we find that 
they contribute about equally to the growth of large 
KBOs, if 

10"^ (9) 



Therefore, growth of large objects will be dominated by 
the accretion of small bodies, if E/cr <^ c?l^. If on the 
other hand E/cr c?l^, then accretion of comparable 
size objects will be the dominant mode of growth pro- 
vided that V remains less than the Hill velocity. Since 
initially Yi/a <^ cc'l'^ it follows that the growth of large 
KBOs was at the beginning dominated by the accretion 
of small bodies, as assumed in equation ([3]). This mode 
of growth continues until E/cr ~ c?!'^ at which stage 
accretion of comparable size bodies starts to compete 
with growth by accretion of small bodies. From then 
on, large KBOs grow in roughly equal amounts by ac- 
creting small bodies and by merging with comparably- 
sized KBOs. As a result E remains constant, since 
its increase due to accretion of small bodies is counter- 
acted by a decrease due to the accretion of large KBOs. 
The radii of large KBOs grows linearly with time (see 
equation [B]). We find it very encouraging that the ob- 
served mass in large KBOs in each logarithmic size in- 
terval estimated from recent Kuiper bel t surveys (e.g. 
Petit et all l200l iTruiillo fc BTott] [2003L iTrujillo et ahl 
2001h is~about 10"'^ of that of the minimum mass so- 
lar nebula (MMSN) (jHavashil [1981 extrapolated to a 
heliocentric distance of 40 AU. The observed mass in 
large KBOs (i? > 50 km) is therefore consistent with 
E/cr ^ 10"'^ and with the hypothesis that they formed 
by coagulation in the Kuiper belt from an MMSN. The 
growth of KBOs ended in this runaway phase with E '-^ 
ot'l'^a ~ 10"'^ cr. In the context of planet formation this 
phase of runaway growth is terminated by the onset of 
oligarchic growth, when each large body dominates the 
stirring in its own feeding zone, allowing E to become 
comparable to a. 

Given this understanding of the KBO growth we now 
show that our choices for the velocity regimes of u and 
V are self-consistent with this picture. Substituting 
Yj/a ^ 10""^ into equation (j4]), we find that u ~ ivn- 
This implies that u is about a few times the Hill velocity 
of the large bodies. We note here that our derivation 
of the expression for u, specifically equating equations 
([21) and ([2]), remains valid even when the large bodies 
contribute significantly to the growth of large KBOs, be- 
cause the growth rates due to the accretion of small and 
large bodies are comparable. Evaluating equation ([6]) we 
find that v ~ Q.Ivh- This implies that the velocity dis- 
persion of large bodies was sub-Hill (i.e. v < vh ) during 
the formation of large 100-km-sized KBOs while that of 
the small bodies was super-Hill (i.e. u > vh)- This con- 
firms that we used the correct velocity regime for u and 
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V in the derivation of the KBO accretion rates above, 
ensuring that our treatment is self-consistent. 
Rewriting equation ([9]) as 



era 



3/4 



constant. (10) 

As we have shown above, S of the largest bodies is con- 
stant with a value of ~ a^^^a. Once large bodies form, 
their mass per logarithmic mass interval is, apart from 
a very brief period (see section \2A^ . preserved. This is 
because bodies smaller than the largest objects, do not 
grow significantly on the growth timescale of the largest 
objects, i.e. the growth is in the runaway regime (see 
section 1^31) . Furthermore, such smaller bodies are also 
not efficiently consumed by larger objects. Their mass 
surface density and size distribution are therefore frozen 
(see section 12. 3p . Thus, S is constant in time, which 
results in 

N{>R)(xR~^, (11) 

since S oc A^(> R)R^. This implies a power-law index 
g = 4 (see equation ([T])) for large KBOs. 

Our work suggests that the growth of large KBOs 
resulted from the accretion of small and large KBOs 
alike and that this mode of growth in the runaway 
regime gave rise to the observed size distribution of 
large KBOs. We confirm this result in section |3] us- 
ing coagulation simulations. Our findings are con- 
sistent with direct observations of the size distribu- 
tion of large KBOs (Bernstein et al. 20041: iFuentes et aD 
120091 : iFraser fc Kavelaarai200&) and also agree with re- 
sults from numerical coagulation simulations that model 
the growth of KBOs carried out by other groups (e.g. 
Eenyon & Luu 1999; Kcnyon 2002). 

2.2. Intermediate sized bodies: Velocity Dispersion and 

Growth 

So far we have only considered two sizes of bodies, 
small ones and large ones. We now turn our attention 
to intermediate size bodies with radii R' , mass surface 
density E' and velocity dispersion v' . 

In sections 12.31 and 12. 4[ we show that S' is of the 
order E. This is because large bodies form by accret- 
ing small bodies until their mass surface density reaches 
E/cr ^ a^^^. From then on, the mass surface density 
remains roughly constant in a given logarithmic mass 
interval. E E' because once larger bodies form, the 
velocity dispersion of the now intermediate sized bod- 
ies quickly grows to super Hill velocities (i.e. v' > vh), 
which implies that such intermediate sized bodies do not 
grow significantly on the growth timescale of the largest 
KBOs and that they are not consumed efficiently by the 
larger objects. The size distribution of such intermediate 
sized bodies and E' therefore remain constant. During 
the brief period over which v' < vh, intermediate sized 
bodies are efBcientl y ac creted by the largest objects, but 
we show in section 12.41 that this phase is very short and 
that it extends over less than a factor of 2 in radius. 
As a result, E' is of the order of, but slightly depleted 
compared to, E. 

2.3. Intermediate sized bodies with v' > vh 

There are three different velocity regimes for interme- 
diate sized bodies that we have to consider separately. 



In the first regime the intermediate sized bodies' velocity 
dispersion exceeds the Hill velocity of the largest bodies, 
i.e. v' > Vh, and the bodies themselves are sufficiently 
small such that their velocity dispersion is not efficiently 
damped by dynamical friction. In other words, the dy- 
namical friction timescale for these intermediate sized 
bodies exceeds the growth timescale of the large bodies. 
The velocity dispersion of these bodies is therefore domi- 
nated by gravitational stirring from the large bodies and 
it grows on the same timescales as the size of the large 
bodies. This implies that v' ^ u. 



V 

VH 



u 

VH 



E 



1/2 



v-1/8 



(12) 



where we substituted Y,/a ^ o?l^ from equation ([9]) in 
the last step. Equating the dynamical friction timescale 
for bodies of size R! to the growth timescale of the large 
bodies and solving for i?', we have 



R a 



vl/4 



(13) 



This implies that v' ^ u for intermediate sized bodies 
with R' ^ a^/'^R ~ O.li? and smaller. For bodies above 
this size, damping by dynamical friction is important. 

In the second velocity regime, which applies for bod- 
ies larger than a^/^R ~ O.li? and that have v' > vh, the 
evolution of the velocity dispersion is dominated by grav- 
itational stirring from the large bodies and damping by 
dynamical friction generated by the small bodies. This 
yields a velocity dispersion that is governed by 

-4 / \ -4 



I dv' ^ E 

il — a 
pR 



dt 



vh 



pR' 



.(14) 

Balancing the stirring and damping rates and substitut- 
ing v'fj = vh{R' /R), we find that v' is given by 



3/4 



-3/4 



a 



1/16 f ^ 
R 



-3/4 



(15) 

where we substituted for u from equation ([4]) and used 
Ti/a '-^ o?!'^ in the last step. Equation (fTSj) yields that 
v' - VH for KBOs with R! - a^l^'^R ~ 0.5 R. Therefore, 
v' > Vh for bodies with R! < 0.5i?. This implies that 
KBOs with radii smaller than ^ 0.5i? have super-Hill 
velocities (i.e., v' > vh) and those with radii larger than 
~ 0.5i? will have sub-Hill velocities (i.e., v' < vh)- 

The growth of intermediate sized bodies that have v' > 
Vh and R' > o^^'^R ^ O.li? is given by 



1 di?' ^ cr /i?'\ / u 
'R'~dt ^ 'pR \r) \^ 



-n 



E' /i?' 



pR \ R J \vh 
(16) 

Substituting for u from equation ([4]) and using again the 
result that E/cr ^ a^^^ we can write the above expression 
as 



R' dt pR\R) pR\R)\vH 

(17) 

where the first term corresponds to growth by the accre- 
tion of small bodies and the second term to growth by 



i?' 



E' /i?' 
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mergers of comparable sized bodies. Since v' > vh, we 
find when comparing the magnitude of the two growth 
terms in equation (jl7p . that the growth is dominated by 
the accretion of small bodies. We assumed that E' is of 
the order of E, which we confirm in section [2.41 Inter- 
mediate sized bodies with i?' < 0.5R therefore grow pre- 
dominantly by accreting small bodies and their growth 
rate is reduced by a factor of {R' /R) compared to bodies 
of size R. This implies that bodies of size R' do not get 
the chance to grow, compared to the growth time scale 
of bodies of size R. 

These intermediate sized bodies do not contribute sig- 
nificantly to the growth of the large bodies. When ex- 
amining the contributions from bodies with v' > vh to 
the growth of bodies of size R, we have 



R at pR pR \vh 



(18) 



where the first and second term correspond to growth by 
merging with bodies of size R and R' , respectively. Since 
the first term in equation ([T5| exceeds the second, large 
bodies grow predominantly by accreting small bodies and 
large bodies. Intermediate sized bodies with v' > vr 
are not important for this growth. Furthermore, such 
intermediate sized bodies are only inefficiently accreted 
by bodies of size R. This is apparent when examining the 
rate of change of their surface density due to accretion 
onto bodies of size R: 



1 dE' 



pR \VH 



(19) 



Comparing equation (|T9)) with equation ([7]) , we find that 
E' does not change significantly on the growth timescale 
of the large bodies, which implies that their surface den- 
sity is not altered due to accretion onto large bodies. As 
a result the mass surface density per logarithmic mass in- 
terval of intermediate sized bodies with v' > vh remains 
constant in time. 

2.4. Intermediate sized bodies with v' < vh 

Intermediate sized bodies that are about half the size of 
the large bodies and larger will have a velocity dispersion 
that is smaller than the large bodies' Hill velocity. The 
evolution of the velocity dispersion of such intermediate 
sized bodies is determined by gravitational stirring from 
large bodies and damping by dynamical friction gener- 
ated by the small bodies. The expression for the evolu- 
tion of v' can be written as 



1 dv' ^ E _2 / 
v' dt pR \vh 



-n—a 

pR' 



.(20) 

Balancing the stirring and damping rates and substitut- 
ing v']^ = vh{R' / R), we can write v' as 



VH 



a ~ \ - 

a 



E\^ /i?'^ 



R 



,1/4 



(21) 



where we substituted for u from equation (|4]) and used 



E/cr 



a' 



3/4 



in the last step. Strictly speaking, we should 



multiply the right hand-side of equation (|211) by a loga- 
rithmic factor, which is given by 3 x Log{R/ R<f), where 



i?* is equal to i?', if v' < v'^. If v' > v'^ then i?* is the 
radius that corresponds to bodies with a velocity dis- 
persion that is less than or equal to v'^. The base of 
the logarithm is the same as the base of the logarithmic 
mass interval over which E is defined. This logarithmic 
factor arises because all bodies larger than or equal to R' 
that have a velocity dispersion less than v'^j contribute 
to the gravitational stirring for q = A. This is a key dif- 
ference compared to the super-Hill velocity regime where 
the stirring is dominated by the large bodies alone. 

The growth of intermediate sized bodies with R' > 
a^/^^R 0.5R is given by 



I dR' ^ cr /R'\ / u 



^-^Facc (22) 



where 



-1 ( v' 



v-3/2 



if v' > v'j_ 
if v' < v'. 



H 



(23) 



H- 



The radius, Rv'^, for which v' ^ v'^j can be found from 
equation (|2ip when including the logarithmic factor and 
it is given by R^'^/R - ai/iS(3 x Log[R/ R^^Jfl^ . The 
v' > v'^ regime given in equation (1221) yields the same ex- 
pression for the accretion rate as given by equation (jl7p . 
The growth in this shown above, is dominated by 

the accretion of small bodies. In the second case, which 
v' < v'jj, the rate of growth of bodies with radii R' is 



1 dR' ^ E fR'\ _3,2 ^ E' 
R' dt pR\R J pR' 



-3/2 



(24) 



where we substituted again for u from equation [4] and 
used Yi/a o?^'^ . The first term corresponds to growth 
by the accretion of small bodies and the second to growth 
by merger of comparable sized bodies. Comparing the 
magnitude of the two growth terms in equation (j24l) 
we find that, unlike in the regime discussed above, the 
growth of bodies with radius R! is dominated by the ac- 
cretion of similar sized bodies rather than by the accre- 
tion of small bodies, again assuming that S' is of the 
order of E, which we show below. We note here that 
similar to the sub-Hill velocity excitation rate, the mass 
accretion rate that corresponds to the accretion of com- 
parable sized bodies, should have been multiplied by a 
logarithmic factor, which is given by 3 x Log\R! jR^'^. 

The maximum value for R' / R^i^ is given by a"^/^^ 2. 
Strictly speaking, we should also have included this log- 
arithmic factor in the derivation of equation ([9]). How- 
ever, since this logarithmic factor only provides a small 
correction it does not warrant to be included here, since 
we have been neglecting factors of order unity through- 
out. Although the overall growth proceeds in the run- 
away regime (i.e, (dLog[R]/dt)/(dLog[R']/dt) = R/R'), 
large bodies with v' < v'^j grow in an orderly fashion 
with respect to each other if q = 4, which implies that 
their radii converge. This leads to a steepening in the 
size distribution at very large radii, but only for bodies 
with radii that have a corresponding velocity dispersion 
such that v' < v'^j, which corresponds to less than a fac- 
tor of 2 in radius. The largest effect this can have is to 
reduce E by a factor of ~ 2, such that E' ~ 0.5E, since 
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orderly growth will cease when S/i? ^ S'/i?'. 

Intermediate sized bodies with v' < vh significantly 
contribute to the growth of large bodies, since 



1 dR 
R~dt 



n—a 

pR 



-3/2 



n — a 
pR 



3/2 



(25) 



where the first and second term correspond to growth 
by merging with bodies of size R and i?', respectively. 
This implies that for q — i.e. E ~ E', large bodies 
of size R grow at comparable rates by accreting bodies 
of their own size and bodies of size R' > 0.5R. Because 
these intermediate sized bodies are efficiently accreted by 
large bodies, their surface density is altered at a rate 



E' dt pR 



(26) 



Comparing equation ((26)) with equation (O , we find that 
E' changes on the order of the growth timescale of the 
large bodies. As a result, the mass surface density per 
logarithmic mass interval of intermediate sized bodies 
with v' < Vh gets depleted, because such bodies are ef- 
ficiently consumed by large bodies. However, since the 
size range is very small, such bodies are only depleted by 
factor of a few before their velocity dispersions become 
super-Hill at which point they are no longer efficiently 
accreted. The orderly growth of bodies with v' < v'^ 
relative to each other and the efficient consumption of 
bodies with v' < vh are responsible for the steepening of 
the KBO size distribution seen at R' > Ry'^ in Figure [H 

3. COAGULATION SIMULATIONS 
3.1. Model 

The aim of our coagulation simulation is to test our 
analytic results outlined above. Its purpose is to cap- 
ture the dominating physical processes that give rise to 
the KBO size distribution. We attempt by no means 
to present the most detailed or precise KBO forma- 
tion simulation, since several su ch works already ex- 



ist in the current literature 



e^ 



IKenvon|[200l IKenvon fc BromlevI 



Kenvon fc (l999l: 



We therefore 

neglect factors of order unity in the accretion, stirring 
and damping rates. In addition, we neglect effects of 
gas damping and possible dynamical stirring from Nep- 
tune. We investigate the KBO growth in a single an- 
nulus located at roughly 40 AU from the Sun with a 
width of about 10 AU. Most of our simulations start 
with a total mass of about 40 Earth masses in small 
planetesimals, which corresponds to tr ~ 0.65 g cm^^. 
This mass surface density is consistent with extrapola- 
tions of the MMSN (Hayashi 1981) to 40 AU, after it has 
been enhanced about 6 fold as required for the forma- 
tion of Uranus and Neptune fe.g. iGoldreich et al.ll2004al : 
iDodson- Robinson fc BodenheimeHl2010D . We follow the 
mass growth and the evolution of the velocity disper- 
sion of the KB Os using Safronov's statistical approach 
(|Safronovlll972f ). Since we are primarily concerned with 
the initial growth phase, we assume here that all colli- 
sions lead to accretion. This assumption is justified dur- 
ing the initial KBO growth and remains justified even 
once the planetesimal's velocity dispersion has been ex- 
cited above its own escape velocity as long as the growth 
timescale of the large KBOs is short compared to the 



planetesimal-planetesimal collision time. This condition 
is generally fulfilled, if the initial planetesimals are of 
the order of a kilometer in size or larger, since, for 1- 
km-sized bodies, the planetesimal-planetesimal collision 
time is comparable to the formation time of Pluto. We 
therefore neglect the effect of destructive planetesimal 
collisions here, but we discuss how they would effect our 
results, had they been important, in section 4. 

Below we give the relevant stirring, damping and ac- 
cretion rates that determine the growth and velocity evo- 
lution in the Kuiper belt. A detailed derivation of the 
these ra tes can, for example, be found in IGoldreich et all 
(|2004bD . The accretion rate between bodies of two dif- 
ferent mass bins is given by 



1 



VHb < '^rel < 



pRb \Ms) 1 _A""b) 

\ a Vrel < VHb 

(27) 

where the subscripts 's' and 'B' correspond to the mass 
bin with the smaller and larger bodies, respectively. The 
number of bodies in a given mass bin is given by N and 
Vrel is the relative velocity, which was approximated by 
Vrel = Max[vsTVB\- Notc, we quote here the sub-Hill ac- 
cretion rate applicable to large objects, since the velocity 
dispersion of the small bodies is super-Hill throughout 
the growth. 

The corresponding velocity evolution for each mass bin 
is dominated by the following processes. The velocity 
dispersion of the mass bin that corresponds to the smaller 
objects is viscously stirred by the object in the larger 
mass bin at a rate given by 



1 dvs „ Es I Q,-2 f iv^ 
Vs dt pRb / X -1 



'^esCB ^ ^rel 

VHb < Vrel < Vesc 



Vrel < VHb 



(28) 

The velocity dispersion of the bodies in the larger of the 
two mass bin is in turn damped by the smaller bodies. 
The damping rate of the large bodies velocity dispersion 
is given by 



1 dvi 



vb dt pRi 



'^escB *^ "^rel 

VHb < Vrel < VescB 
Vrel < VHb 

(29) 

We implemented the above equations in our coagulation 
code and followed the growth and velocity evolution of 
the various mass bins. 



3.2. Results 



3.2.1. 



KBO Growth and Velocity Evolution 

As initial conditions we start with an equivalent of 
about 40 Earth masses, in 1-km sized planetesimals with 
an initial velocity dispersion of 3 times their Hill veloc- 
ity. We follow the KBO growth and velocity evolution of 
the different mass bins for 70 Myrs. The results of our 
coagulation simulation are shown in Figures [I] [21 and [31 

Figure [1] shows the evolution of the cumulative mass 
distribution, the cumulative number distribution and the 
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velocity dispersion during KBO growth as a function 
of time. It is apparent from Figure [T] that the growth 
and velocity evolution become self-similar once S ~ 
o?l'^(T ~ 10~^cr, i.e. the shape of the size-distribution 
and velocity-distribution remains unchanged, while the 
maximum KBO size continues to grow. 

The two upper panels of Figure [T] and Figure [2] show 
that the size-distribution of large KBOs (R > 50 km) in- 
deed follows a power-law with g = 4 as predicted by our 
analytic treatment in section [2] This implies a roughly 
equal amount of mass per logarithmic mass bin for large 
KBOs. Moreover, the middle panel of Figure [T] and Fig- 
ure[2]show that the mass ratio of large to small KBOs, i.e. 
S/fj, found in our coagulation simulation agrees very well 
with our analytic prediction that S/cr ^ a^'^/'* ^ 10"'^. 
The simulations confirm our analytic results and sug- 
gest that the total mass in large objects that we see in 
the Kuiper belt is not arbitrary but an outcome of the 
KBO growth and that it is roughly 10~^ of the initial 
planetesimals mass. This result is in excellent agree- 
ment with the actual observed mass in large KBOs and 
formation from a MMSN type disk. Our work there- 
fore suggests that the Kuiper belt did not contain two 
to three order of magnitude more mass in large KBOs 
as has be en proposed by some models (jWeidenschillind 
[2002; Tsig anis et al|[2005l ). 

Figure [3] shows a comparison of the velocity dispersion 
from our coagulation simulation and our analytic results 
derived in section [2j It displays very good agreement 
between the velocity dispersion that we derived analyt- 
ically for the various size regimes, and the results from 
our coagulation simulation. 

Figures [l][3] show that our analytic work captures the 
essential features of KBO growth and that analytic the- 
ory and the numerical coagulation results are in excellent 
agreement. We are able to successfully explain the slope 
{q = 4) and amplitude (E ~ a~^/'^a) of the large KBO 
size distribution and the evolution of the velocity disper- 
sion in the various velocity and size regimes. 

3.2.2. Initial Planetesimal size distribution and sizes 

We performed an additional set of coagulation simula- 
tions with different initial conditions. In the first set we 
started with most of the mass in small 1-km-sized ob- 
jects, just as before, but added a second population of 
larger 10-km-sized KBOs that contained 10^'' of the total 
mass. All bodies were started with a velocity dispersion 
equal to 3 times their Hill velocity and we followed their 
growth and velocity evolution. Figure |3] shows the result 
of this growth (points) and a comparison with the cumu- 
lative mass distribution at various times for KBOs that 
grew from a single population of 1-kni sized planetesi- 
mals (lines). The similarity between the two distribu- 
tions is striking. The same power-law for the large KBO 
size distribution emerges and the mass ratio in large and 
small KBOs becomes also the same in both simulations. 
These results highlight the power of the 'two groups ap- 
proximation' that we used to derive the analytic results 
in section [2] and validate our assertion that the growth 
of large KBOs develops towards a state where small and 
large bodies contribute about equally to the growth. The 
overall growth timescale to reach Pluto-mass objects dif- 
fers in these two simulations. Starting from 40 Earth 




Radius [km] 




1 10 100 1000 

Radius [km] 

Figure 1. Evolution of the cumulative mass distribution, the cu- 
mulative number distribution and the velocity dispersion during 
KBO growth as a function of time at 3 X 10' years (dashed or- 
ange line), 4 X 10^ years (dotted green line), 5 X 10^ years (dot- 
dashed turquoise line), 6 X 10^ years (long-dashed blue line), and 
7 X 10^ years (solid purple line). The growth of the large KBOs 
{R > 50 km) becomes self-similar from 4 X 10^ years onwards, i.e. 
the shapes of the size-distribution and of the velocity-distribution 
remain unchanged but the KBO size and velocity continue to grow 
with time. The slope of the large KBO size distribution is g ~ 4. 
A power-law index of q = 4 corresponds to a horizontal line in the 
top panel of this figure. The thin black line in the middle panel 
of this figure represents the KBO size distribution with q = 4, 
as predicted by our analytic theory and its expression is given in 
the top righthand corner. The bottom panel shows the evolution 
of velocity dispersion during KBO growth. The escape velocity 
and Hill velocity (assuming a KBO material density of 1 g/cm^) 
are given as a function of size by the upper and lower thin black 
line, respectively. We refer the reader to Figure |3] for a detailed 
comparison between our analytic theory and the velocity evolution 
found from our simulations. Our analytic theory and the numerical 
coagulation results are in excellent agreement. 
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Radius [km] 



Figure 2. KBO size distribution as a function of radius at 
5 X lO'jears from the same coagulation simulations as shown in 
Figure m The y-axis corresponds to the mass in a given Log2 
mass interval. The total mass in small planetesimals is given by 
a and its values is given by the orange dashed line drawn below 
its symbol. Similarly, the mass in large KBOs in a given Log2 
mass interval is denoted by S and its predicted value from section 
[2] is given by orange dashed lines drawn below it. The value for S 
from our analytic work and numerical coagulation simulation agree 
within a factor of 2. This agreement can even be improved, if we 
account for the fact that all logarithmic mass bins with v' < vn 
contribute to the growth of th e la rgest bodies about equally (see 
discussion following equation I I24I1 in section O , which would re- 
duce the predicted value of S somewhat. The dashed vertical lines 
mark the radii separating the different velocity dispersion regimes, 
as discussed in section [2] 



1000 ■ 
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Radius [km] 

Figure 3. Velocity dispersion as a function of radius at 5 X 
10 years. The results from our coagulation simulation are given 
by the black points, the dashed orange lines are our analytic pre- 
dictions for the velocity dispersion with their equations given above 
each segment. The dashed vertical lines mark the radii separating 
the different velocity dispersion regimes, as discussed in section |2| 



masses in small planetesimals, it takes in our simulations 
about 70-80 Myrs to form Pluto-sized objects if initially 
all KBOs are 1-km in radius but only 40-50 Myrs if, in 
addition to the 1-km-sized objects, there also existed a 
small population of 10-km-sized KBOs. The difference 
in the growth time results from the timescale it takes 
1-km-sized objects to grow into 10-km sized bodies, if 



none are already present. When only comparing the 
growth timescale from 50-km-sized objects to 1000-km 
sized objects, we find that they are the same in the two 
scenarios (see figure U]) , confirming that it is the initial 
growth of objects to tens of kilometers in size that give 
rise to the overall difference in Pluto-formation time be- 
tween the two simulations. The final shape of the KBO 
size-distribution and the total mass in large KBOs (as a 
fraction of the initial mass) is therefore independent of 
the initial size distribution of the planetesimals. We also 
find the same KBO size distribution and overall mass in 
large KBOs, if the initial planetesimals have a power-law 
size distribution with various power-law indexes and radii 
ranging from 1 to 10-km. In addition, we confirm that 
our results are independent of the details of the initial 
planetesimal velocity dispersion, as long as it is below 
their escape velocity. 




10 100 
Radius [km] 

Figure 4. Direct comparison of KBO growth from two different 
initial conditions. The lines correspond to starting with only 1-km- 
sized bodies. The points represent the KBO growth that resulted 
when starting with the same total mass as before, but in addi- 
tion to having most of the mass in 1-km-sized objects we added a 
second population of 10-km-sized KBOs that contain 10""^ of the 
total mass. The cumulative mass distribution of the small objects 
is omitted here, since it is the size-distribution of large KBOs that 
wo want to compare here. The agreement between the two simula- 
tions is striking. The shape, amplitude (i.e. the ratio of S/cr) and 
power-law slope is therfore not a result of the initial conditions but 
represents an 'equilibrium state' that the system evolves to. The 
dotted green line corresponds to a time of about 3 X lO'^ years since 
the start of the simulation, whereas the green points correspond to 
only about 10^ years since the start of the simulation. Due to the 
different initial conditions it took a different amount of time to 
grow to 50-km-sized KBOs, but the growth timescale once these 
50-km-sized KBOs have formed (green dashed line and points) be- 
comes the same. In other words, it takes the same amount of time 
to grow from 50-km-sized KBOs to Pluto-sized objects irrespective 
of the initial conditions, once 50-km-sized KBOs have formed. The 
blue dashed line and points and the red solid line and points corre- 
spond to 1.5 X 10^ years and to 5 X 10^ years, since the formation of 
objects that are represented by the green dashed lines and points, 
respectively. 



We also performed a second set of coagulation simula- 
tions, which we started with the same total planetesimal 
mass but with an initial planetesimal size of 500 m in 
radius rather than 1 km. Figure |S] shows that the shape 
of the size distribution does not depend on the initial 
planetesimal size. This result remains valid as long as 
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the initial planetesimal size is large enough such that the 
planetesimal collision timescale exceeds the KBO growth 
timescale. 




10 

Radius [km] 



100 



Figure 5. Comparison of KBO growth from two different initial 
planetesimal sizes. In each simulation, all planetesimals started out 
with one size. The solid and dashed line correspond to an initial 
planetesimal size of 500 m and 1 km, respectively. The mass distri- 
bution is shown after a growth time of ^ 5 X 10' years. The y-axis 
corresponds to the mass in a given Log2 mass interval. The shape 
of the size distribution does not depend on the initial planetesimal 
size, as long as the planetesimal-planetesimal collision timescale 
exceeds the KBO growth timescale, such that dynamical cooling 
and possible mass loss due to frequent planetesimal collisions can 
be neglected. 



These resu lts confirm earlier findings by 

iKenyon fc Luul (|1999( ). who note that the power- 
law slope for large objects is remarkably independent 
of the input parameters and initial conditions. Our 
work offers an explanation of why this is so. Since the 
size distribution of large KBOs reaches an equilibrium 
state that evolves self-similarly, the signature of the 
initial planetesimal size-distribution and their initial 
velocity dispersion are erased. The slope of the KBO 
size distribution and the formation timescales that we 
find in our simulations for large KBOs is in agreement 
with results from previous coagulation simulations 
(iKenvon fc Luulll999l : Ikenvonll2002t iKenvon fc BromlevI 

3.2.3. Different values of a 

Since our analytic result for the mass ratio between 
large and small KBOs (i.e., S/ct ~ a^!^) can be ex- 
pressed as a function of a alone, we performed an addi- 
tional set of coagulation simulations for a value of a that 
was decreased by a factor of 100. Figure [6] shows the 
coagulation results for a = 10"'' (dashed lines), which 
is roughly the value for the Kuiper belt, and a — 10~^ 
(solid lines). We increased the initial planetesimal mass 
in the a — 10~^ simulation by a factor of 10^-^ to speed 
up the planetesimal growth. Due to this mass increase, 
the resulting mass per logarithmic mass interval in large 
bodies should be the same in both simulations. This 
is indeed what we see in Figure [SI confirming our an- 
alytic expression for E. Figure [6] shows also that the 
equilibrium growth state, where growth by accretion of 



small bodies is comparable to the growth due to similar 
sized mergers, is reached later for a = 10^^ compared 
to a = 10"**. We confirmed that this is not caused by 
the increased planetesimal mass that was used in the 
a = 10^^ simulation. 
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Figure 6. Comparison of our coagulation results for two different 
values of a. The dashed lines correspond to ai = 10""*, which is 
roughly the value for the Kuiper belt, and the solid lines correspond 
to a2 = 10~®, which could, for example, represent a debris disk 
around a solar-mass star with a semi-major axis about 100 times 
larger than the Kuiper belt. We increased the initial planetesimal 
mass in the «2 = 10~® simulation by a factor of 10^'^ to speed up 
the planetesimal growth. The y-axis corresponds to the mass in a 
given Log2 mass interval; ci and (T2 are the sum of the total mass 
in small planetesimals. Si and S2 correspond to the mass in large 
KBOs in a given Log2 mass interval and their predicted values 
from section^are given by the lines drawn below their respective 
symbols. The agreement between theory and simulation for both 
systems is clearly shown in Figure |6] 



4. COLLISIONS AND OTHER CONSIDERATIONS 

We have so far neglected collisions among the small 
bodies. In this section, we first derive the planetesimal- 
planetesimal collision timescale and show that plan- 
etesimal collisions can be neglected during the KBO 
growth for planetesimal sizes of about a kilometer and 
larger. We then discuss how our results would be al- 
tered, if planetesimal collisions had been important dur- 
ing the KBO growth. In addition, we examine the ef- 
fect of semi-collisional a ccretion (Sar i fc Goldreicbll2006l : 
ISchlichting fc Sarill2007[ ) and frequent binary mergers on 
the growth of KBOs. 

The planetesimal-planetesimal collision time for ob- 
jects with radius, r, and with a velocity dispersion that 
has been excited above their escape velocity is given by 



TColl 



n- 



(30) 



Comparing equation (l30l) to the Pluto formation 
timescale, which is given by the inverse of equation ([S]) 
we find that 



TColl 



.-3/2 



which is 



Tpiuto XCrRphito 

1 for 1-km-sized planetesimals and E 



(31) 



1-km is about the size that planetesimals are 
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expected to have, if they formed by gravitational insta- 
bility without dissipation of internal angular mom entum 
(|Goldreich fc WardlHOTl IGoldreich et alJl2004bn . This 
implies that if KBOs grew from a population of 1-km 
sized planetesimals then mutual planetesimal colhsions 
become important only once KBOs comparable to the 
size of Pluto have formed. 

If planetesimal-planetesimal collisions are important, 
then they would, most likely, break-up the planetesimals 
and lead to a damping of the small bodies velocity dis- 
persion, which would therefore no longer grow as R, as 
we assumed in section [21 Instead, it would be set by 
balancing the dynamical stirring from large bodies by 
collisional cooling. This yields 



u 

VH 



(32) 



for u < vh, where s is the radius of the small bodies 
(IGoldreich et al.l l2004bl ) . We note here that u < vh is 
most likely the relevant case to consider, since u is only 
slightly super-Hill when collisional cooling is neglected 
(see section[2]). In this case, the growth of KBOs is given 

by 



IdR _ifu\' ^E 

R dt pR \VH J pR 



3/2 



(33) 



where we assumed that a^^^vn < u < vjj- Substituting 
for u from equation p2[) and comparing the two growth 
rates in equation (j33|) , we find that the two growth rates 
become equal when 



E 
a 



— a 



1/2 



10 



_i /i?/1000 km 
\ s/1 cm 



1/2 



(34) 



This implies, that in this case, the growth of large bodies 
is dominated by the accretion of small bodies until more 
mass is converted into S, compared to the case without 
frequent planetesimal collisions. The total mass in large 
KBOs and the slope of their size distribution are there- 
fore likely to be different from our results in section [2] 
and [21 if collisions were important. The KBO growth 
would be halted, if collisions lead to the onset of a colli- 
sional cascade and if the small collisional fragments are 
efficiently lost from the Kuiper belt. 

So far, we have treated the planetesimal accretion as 
collisionless, meaning that collisions among planetesi- 
mals can be neglected while they are inside the Hill 
sphere of a growing KBO. This assumption is valid, if 
KBOs formed from km-sized planetesimals. However, 
if planetesimals are of the order of a meter in size or 
smaller, either because they formed small, or because 
they were broken into small pieces, then they are more 
likely to collide with each other inside a KBO's Hill 
sphere than to accrete o nto the growing KBO directly 
(ISchlichting fc Saril [20071) . We call this semi-collisional 
accretion. Such planetesimal collisions inside the Hill 
sphere lead to the formation of an accretion disk around 
the KBOs. If KBOs grew by semi-collisional accretion 
than their growth could have been very fast, because the 
effective radius for accretion, in this case, is of the order 
of their Hill sphere. The KBO growth would be domi- 
nated by the accretion of small planetesimals for much 
longer in the semi-collisional regime, compared to the 



growth scenario discussed and investigated in sections [5] 
and[3l 

Finally, the growth of KBOs might also have been 
aided by merging comparable mass Kuiper belt bina- 
ries. As we discuss in section [5| in detail, a significant 
fraction of KBOs reside in comparable mass binary sys- 
tems and such systems likely formed by dynamical pro- 
cesses such as by three body capture, or by bin ary forma- 
tion aided by dynamical friction ( Wcidensch illinel 120021: 
IGoldreich et al.''200 ^lFunato et al.ll2004tlAstakhov et al.l 
2005; Lcc ct al. 2003). Once formed, the mutual semi- 
major axis of such binary systems starts to shrink due to 
dynamical friction gen erated by the small planetesimals 
(jGoldreich et al.ll2002[ ). This eventually leads to merging 
between the binary components. This channel of growth 
could be important, because the binary formation time 
and the time it take for a binary to spiral in until contact 
occurs are c omparable to the tim escale on which large 
KBOs grow (IGoldreich et al.ll200l) . If merging of binary 
components is important in the overall KBO growth then 
the growth rate would be enhanced compared to equation 
([7]). This results in a reduced mass surface density for 
large KBOs, but the power-law slope of the large KBO 
size distribution will remain unchanged. 

5. IMPLICATIONS FOR KUIPER BELT BINARIES 

Our coagulation results also have interesting implica- 
tions for the formation of Kuiper belt binaries. A sig- 
nificant fraction of KBOs are part of a binary system. 
The binary fraction varies for different dynam ical classes 
and it is highest in the cold classical belt (jNoll et al.l 
I2008cf ). where it is about 30%. More than 70 bina- 
ries have been discovered in the Kuiper belt to date 
and their number continues to rise. High-mass ratio 
binary systems, including Pluto/Charon, likely formed 
in a collision and subsequent tidal evolution. However, 
the majority of Kuiper belt binaries consist of compa- 
rable mass companions with wide separations and have 
too much angular momenta to have formed by the same 
mechanism. Instead these systems most likely have a 
dynamical origin like, for example, binary formation by 
three body capture, or binary formation aided by dy- 
namical friction ([W eidenschilling 2002; Goldrcich et al 



2002HFunato et aLll2004tiAstakhov et al.M2005:,Lee et al 



2007f l. The Hill sphere, the region interior to an object's 



Hill radius, sets the maximum phase space available for 
binary formation by such dynamical processes. Dynam- 
ical binary formation scenarios take advantage of the in- 
creased size of the Hill radius, which is more than an 
order of magnitude larger in the Kuiper belt than for 
similar sized objects in the Asteroid belt and therefore 
make the Kuiper belt the ideal place for the formation 
of wide, comparable mass binaries. 

Dynamical binary formation scenarios require v' < 
Vh for efficient formation, since binary formation rates 
quickly exceed the age of the solar system, once the ve- 
locity dispersion significantly exceeds the Hill velocity 
(INoll et al.l 120083 ISchlichting fc Sarill2008l) . Since KBO 
growth and binary formation occurs concurrently, we can 
predict the maximum mass ratio for such binary sys- 
tems because we know the range of KBO radii that have 
sub-Hill velocities at any given stage during the growth. 
From equation (|15p we have v' < vh for KBOs with 
radii greater than R' ~ a^^^'^R 0.5 R. This im- 
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plies a maximum mass ratio of 



^1/4 



10 between 



the primary and secondary and translates into a max- 
imum magnitude difference between the binary compo- 



nents, assuming similar albedo, of about A,; 



1.7. 



This provides an explanation for the strong clustering 
in Amag of binaries in the cold cl assical bel t. This 
clustering was first pointed out by iNoU et al.l (|2008bf ) 
and is shown in Figure [T] Figure [7] shows that all bi- 
naries in the cold classical belt have A^ag < 1.7, i.e. 
they lie in the region enclosed by the solid red lines. 
Our coagulation work therefore suggests that the ob- 
served binaries in the classical belt formed by dynam- 
ical processes during the growth of the KBOs them- 
selves and that this formation most likely took place in 
situ. The cold classical belt is therefore not only distinct 
in its inclination d i stribu tion (|Brown 2001; Elliot ct al] 
12005^, 'Gulbis ct al.' '2010), co lor (iTe glcr & Romanishhl 
(2000; Trujillo & Brown 2003 iPet it et al. 200i) and bi- 
nary fraction (lStcphcns.^^NoiJl2006) , but also in terms 
of the type of its binaries, i.e. they are exclus ively com- 
parabl e mass ratio binaries with Amag < 1.7 l)Noll et al.l 
l2008al |bl). Figure [7] shows that, in contrast to the cold 
classical population, other dynamical classes have bina- 
ries with a wide range of mass ratios and they there- 
fore most likely contain dynamically and collisionally 
formed binary systems. The fact that the other dynami- 
cal classes have comparable mass binaries as well as sys- 
tems with a wider range of component sizes is consistent 
with the idea that they are a superposition of two original 
populations, namely a dynamical cold (low inclination) 
population that formed close to its current location and a 
dynamical hot (high inclination) population that formed 
closer to the Sun. Kuiper belt binaries may therefore 
prove to be useful probes for entangling the two original 
dynamical populations, if they existed (Murray-Clay & 
Schlichting 2010, submitted). 

6. DISCUSSION & CONCLUSIONS 

We carried out an analytic and numerical investiga- 
tion of runaway growth of protoplanets with a focus on 
the Kuiper belt. Since, the Kuiper belt is a remnant 
of the primordial Solar system where planet formation 
never reached completion, it contains some of the least 
processed bodies. As a result, it provides a snapshot of 
earlier stages of planet formation, which were erased else- 
where in the Solar system where planet formation pro- 
ceeded to completion. Our results for runaway growth, 
which we summarize below, therefore, do not only apply 
to the Kuiper belt but also to planet formation before 
the onset of oligarchic growth. In addition, our findings 
apply to proto planetary growth in debris disks around 
other stars, as long as the protoplanets' growth time 
is shorter than the planetesimal-planetesimal collision 
time, such that dynamical cooling and possible mass loss 
due to frequent planetesimal collisions can be neglected. 

In this paper, we presented analytic work that de- 
scribes the growth of KBOs, the evolution of their ve- 
locity dispersion, and that provides insights into the un- 
derlying physical processes that give rise to the KBO 
size distribution. Our work successfully explains the ob- 
served slope of the KBO size distribution as well as the 
total mass that is present in large KBOs today. In addi- 
tion it predicts the maximum mass-ratio of Kuiper belt 
binaries that formed by dynamical processes, which ex- 




Figure 7. The Hv magnitudes of both components of Kuiper 
belt binaries, labeled by their respective dynamical class according 
to the classification by Gladman et al. (200^). The black dotted 
lines, correspond to increments of 1 magnitude difference between 
the binary components. The region enclosed by the solid red lines 
ranges from equal magnitude binary components to a maximum 



magnitude difference of A„ 



1.7. This corresponds to the 



maximum size difference that Kuiper belt binary components with 
similar albedos can have as predicted by this work, if binary for- 
mation preceded by dynamical processes (see section[5]for details). 
All binaries in the cold classical belt lie in the region enclosed by 
the solid red lines. In contrast to the cold classical population, 
the binary components of other dynamical classes span a wider 
range in mass ratios. The plotted binary data was obtained from 
the follo wing references: Noll ct al. (2008a) and references therein, 
INoU et a l. (2008c) and Lin et al. (2010). We calculated the Hy 
magnitude of each binary component from the magnitude differ- 
ence between the primary and secondary and from their combined 
Hv magnitude as given in the references above. Pluto and Eris, 
which are both not part of the classical belt, are not shown on 
this plot, since they have Hy magnitudes brighter than and fall 
therefor to the left of this figure. 

plains the observed clustering in binary companion sizes 
that is seen in the cold classic belt. We confirmed our 
analytic results with numerical coagulation simulations. 

We find that the KBO growth proceeds as follows. Ini- 
tially all the mass resides in small planetesimals and 
large KBOs start to form by accreting small planetes- 
imals. This growth continues until growth by merg- 
ing with comparable sized KBOs become comparable to 
growth by accreting small bodies. We show that this 
condition sets in when S/cr ^ a'^/'* 10""^. From that 
time onwards, the growth and the evolution of the veloc- 
ity dispersion become self-similar and E remains roughly 
constant, since the increase in S by the accretion of small 
planetesimals is balanced by a decease due to the accre- 
tion of large bodies. We showed that this mode of growth 
leads to a KBO size distribution with power-law index 
(7 = 4. This is in good agreement with observations of the 
Kuiper belt size distribution, which is well described by 
a power-law with an in dex that is consistent with q = 4 
within la (e.g. [T rujillo et al.ll20Qll:lBernsteiii et al.l l2004l: 
'Fuent es fc Holmanl l200'8l : iFraser et al.l 120081 ). A single 
albedo is assumed for all sizes when converting from the 
observed magnitude distribution to the KBO size distri- 
bution. Therefore, possible albedo variations as a func- 
tion of size could introduce a significant uncertainty in 
the estimate of q. The best fit value for the power-law 
index is typically found to somewhat exceed q = A. We 
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note that this could be due to fitting the high end of the 
size distribution that corresponds to the largest KBOs, 
which is somewhat steeper than the rest of the distri- 
bution (see figure [T]). This steeper end of the KBO size 
distribution is due to the fact that the largest objects 
did not have enough time to grow to their steady-state 
abundance. 

If KBOs formed by coagulation from km-sized plan- 
etesimals, then there could not have been significantly 
more mass in large KBOs than what is observed to- 
day, unless the MMSN was initially enhanced by sev- 
eral orders of magnitude. This result is in good agree- 
ment with the current mass in large KBOs, since it is 
about E/cr - a3/4 _ MMSN that was en- 

hanced by a factor of a few, as required for the for- 
mation of Uranus and Neptune. The observed mass in 
large KBOs and their size distribution therefore support 
the hypothesis that KBOs formed from a MMSN type 
disk by coagulation from kilometer sized planetesimals. 
The growth of KBOs ended in this runaway phase with 
E ^ a'^/'^a ^ 10~'^CT. In the context of planet formation 
this phase of runaway growth was terminated by the on- 
set of oligarchic growth when each large body dominates 
the stirring in its own feeding zone. During oligarchic 
growth E continues to grow until it becomes comparable 
to a. 

Our understanding of the growth, the value of E/cr and 
the size distribution of large bodies differs from previous 
works that investigated the details of runaway growth 
(iMakino et all 119981: iLei l2000t iMalvshkin fc GoodmanI 
120011) . This is mainly due to the fact that these works 
did not take into account the simultaneous evolution 
of the velocity dispersion during the growth but as- 
sume d either equipartiti o n between the different mass 
bins ()Makino et al'l 119981: IMalvshkin &: GoodmanI I2001D 
and / or a constant velocity dispersion for all sizes 

Hoop). 

Since efficient binary formation by dynamical processes 
can only proceed among KB Os tha t have sub-Hill ve- 
locity dispersions tNoJl ct aU l2008al: ISchlichting fcSari l 
120081) and because we follow the evolution of the KBO 
velocity dispersion as well as their growth, we can predict 
the maximum mass ratio for such binary systems. Our 
work yields a maximum binary mass ratio of a~^^^ ~ 10, 
or, in other words, a maximum magnitude difference be- 
tween binary components, assuming similar albedos, of 
A 1.7. This explains the clustering in A^ag of all 

the observed Kuiper belt binaries that are part of the 
cold classic belt, which make up about half of all known 
binary systems in the Kuiper belt. 

There exists some tentative observational evidence 
that the power-law size distribution might be different 
for KB Os with inclinations greater and smaller than 
^ 5° (iBer nstein et all 12004 iFuentes fc Holman' '2008| 
iFraser et a l. 2010). Howe ver, Eascr & Kavclaars (200S>i) 
and Fucnt es et al.l (2010) found no conclusive evidence 
supporting such a difference. It will be interesting to 
see the results of future KBO surveys that address this 
question. Should there indeed exist a difference between 
various dynamical classes then this could have very in- 
teresting implications for the formation and collisional 
evolution of the Kuiper belt. 
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